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We report the direct observation of the excited L — 1 state B*2 in fully reconstructed decays to 
B'^K~ . The mass of the 5^2 meson is measured to be 5839. 6± 1.1 (stat.) ± 0.7 (syst.) MeV/c^, and 



its production rate relative to the B^ meson is measured to be [1.15 ±0.23 (stat.) ± 0.13 (syst.)]%. 
PACS numbers: 12.15.Ff, 13.20.He, 14.40.Nd 



To date, the detailed spectroscopy of mesons contain- 
ing a b quark has not been fully established. Only the 
ground J^ = 0" states B+, B", B°, B+ and the excited 
1~ state B* are established according to the PDG [lil. 
Previous studies of excited (bs) states have been carried 
out using inclusive final states with limited precision and 
ambiguous interpretation ^. The properties of (bs) ex- 
cited states, and comparison with the properties of the 
(bu) and (bd) systems, provide tests of various models of 
quark bound states and are important for their continu- 
ing development. 

This Letter presents the observation of the pro- 
cess B*2 -^ B^K~ with exclusively reconstructed _B+ 
mesons, using a data sample corresponding to 1.3 fb^^ 
integrated luminosity collected with the DO detector at 
the Fermilab Tevatron collider during 2002-2006. Charge 
conjugated states are implied throughout this Letter. 

Quark models predict the existence of four P-wave 
(L = 1) states in the (bs) system: two broad resonances 
iB*Q and i?*^) and two narrow resonances (-B^i and B*2) 
[^. The broad resonances decay via S-wave processes 
and therefore are expected to have widths of a few hun- 
dred MeV/c^. Such states are difficult to distinguish, 
in effective mass spectra, from the combinatorial back- 
ground. The narrow resonances decay via Z?-wave pro- 
cesses {L — 2) and should have widths of approximately 
1 MeV/c^ ^4]. If the mass of the B^j (J = 1, 2) is large 
enough, the main decay channel should be B^j -^ B'-*^K, 
since the -B^tt channel is forbidden by isospin conserva- 
tion. A recent result by the CDF collaboration reports 
the observation of two narrow resonances consistent with 
the Bsi and B*2 states [5]. 

The search for narrow Bsj mesons is performed by 
examining events with B^^*^ K~ decays. This sample 
includes the following decays: 



» 

'si - 


^ B*+K- 


, B*+ - 


-^B+r, 


(1) 


1* 

s2 " 


^ B*+K- 


,B*+- 


-^B+r, 


(2) 


's2 


-> B+K-. 






(3) 



The direct decay Bgi -^ B^K^ is forbidden by conser- 
vation of parity and angular momentum. In decays ([1]) 
and ^, the photons from the i?*+ decay have energy 
E{-i) = (45.78 ± 0.35) MeV 1]. These photons are not 
reconstructed in this analysis, so that for such events 
the invariant mass of the reconstructed decay products 
is shifted down by E(^). 

The DO detector is described in detail elsewhere [6|. 
The detector components most important for this anal- 
ysis are the central tracking and muon systems. The DO 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 



located within a 2 T superconducting solenoidal mag- 
net, with designs optimized for tracking and vertexing 
at pseudorapidities |?7| < 3 and |?7| < 2.5, respectively 
(where rj = — ln[tan(6'/2)], and 6 is the polar angle mea- 
sured with respect to the proton beam direction). The 
muon system is located outside the calorimeters and has 
pseudorapidity coverage JTyj < 2. It consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar layers 
behind the toroids jlj. 

The data for this analysis were selected without any 
explicit trigger requirement, although most events sat- 
isfy inclusive single-muon triggers. The B^ mesons are 
reconstructed in the exclusive decay B'^ -^ J/ipK~^ with 
J/'0 decaying to fi^ fi^ . The selection procedure used is 
exactly as described in Ref. [9| . All B mesons with mass 
5.19 < M{B+) < 5.36 GeV/c^ are used, which yields a 
sample of 20915±293 (stat.) ±200 (syst.) B+ candidates. 

For each reconstructed B~^ meson, an additional track 
with transverse momentum (Pt) above 0.6 GeV/c and 
charge opposite to that of the _B+ meson is selected. This 
track is assigned the kaon mass. 

For any track i, the significance Si is defined as Si = 
y[e77^^£70F~H^L7^^^L)P, where er (el) is the projec- 
tion of the track impact parameter on the plane perpen- 
dicular to the beam direction (along the beam direction) , 
and cr(e7') [cr(eL)] is its uncertainty. Since the B^j mesons 
decay at the production point, the additional track is re- 
quired to originate from the primary vertex by applying 
the condition on its significance Sk < v6. The primary 
vertex is defined using the method described in Ref. [8(. 

For each combination satisfying the above criteria, the 
mass difference AM = M{B+K-) - M{B+) - M{K-) 
is computed. The resulting distribution of AM is shown 
in Fig. [11 

Of the three decays (UHl]) through which the Bgj states 
can reach the ground state -B"*", one or more may be kine- 
matically suppressed if the excited state mass is smaller 
than the mass of the decay products. From inspection 
of Fig. [H there is a single region of excess events above 
the background at AM = 67 MeV/c^, therefore the fit is 
based on the hypothesis that only one decay channel is 
observed. From kinematic considerations it follows that 
this is the highest energy transition, i.e. i?*2 -^ B^K~ . 
Alternative hypotheses are discussed later. 

Since the decay i?*2 -^ B^K~ occurs very close to 
the threshold AM = MeV/c^, its width F should be 
around 1 MeV/c^ 4^ . Because this is much less than the 
detector resolution, which is of order 6 MeV/c^, the fit is 
insensitive to values of F below 6 MeV/c^, and F is fixed 
at 1.0 MeV/c^. A systematic uncertainty is assigned to 




M(B*K") - M(BO - M(K") (MeV/c^) 

FIG. 1: Invariant mass difference AM = M{B+K~) - 
M{B'^) — M{K~) for exclusive B decays. The line shows 
the fit described in the text, with signal and background con- 
tributions also plotted separately. 



this choice of F by fitting with a selection of small widths 
in the range to 2 MeV/c^. 

Based on the above, the experimental distribution un- 
dergoes a binned maximum likelihood fit to the following 
function 



F(AM) = i^,ig(AAf) + i^bckg(AM), 
F3ig(AM) = N-D{AM,Ao,r) 



(4) 



In these equations, Aq is the central position of the 
resonance, i.e. MiB*^) - M{B+) - M{K-), F is the 
5*2 width, and N gives the total number of observed 



Bt 



B^K decays. The background is parameterized 



by a modified power-law function: 

Fbckg(AM) = c{AM)''+d{AM), 



(5) 



where the parameters c, d and k participate in all fits. 

The function D(AM, Ao,F) in Eq. (I3|) is the convolu- 
tion of a relativistic Breit-Wigner function with the ex- 
perimental Gaussian resolution in AM. The width of res- 
onances in the Breit-Wigner function takes into account 
threshold effects using the standard expression [l|, [lO| for 
L = 2 decay. 

The detector resolution function is determined from 
Monte Carlo simulation. All processes involving B 



mesons are simulated u sing the EvtGen generator [11 1 
interfaced with PYTHIA I4I , followed by full modeling of 
the detector response with GEANT [13| and event recon- 
struction as in data. The difference between the recon- 
structed and generated values of AM is parameterized 
by a double-Gaussian function, with the width ai ((72) 
of the narrow (wide) Gaussian set to 2.7 MeV/c^ (6.2 
MeV/c^), and the normalisation of the narrow Gaussian 
set to 1.2 times that of the wide Gaussian. Studies of 
the B+ ^ J/ipK+ and D*+ -^ D°tt+ decays show that 
simulation underestimates the mass resolution in data by 



wlO%. Therefore, the widths of the Gaussians which pa- 
rameterise the Bgj resolution are increased by 10% to 
match the data, and a 100% systematic uncertainty is 
assigned to this correction. 

Using a fitting range of < AM < 150 MeV/c^, cov- 
ering 50 bins, a binned maximum likelihood fit is per- 
formed. The following parameters of B*2 are obtained: 

Ao = M{Bl2) - M{B+) ~ M{K~) 

= 66.7 ±1.1 (stat.) MeV/c^, 
N = 125 ± 25 (stat.) events. (6) 

Without the BI2 signal contribution, the log likelihood of 
the fit decreases by 13.4, which implies that the signal is 
observed with a statistical significance of more than ba. 
To convert the Aq result into a mass measurement on 
5*2, the PDG values of the B+ (5279.1 ± 0.5 MeV/c^) 
and R- (493. 677± 0.013 MeV/c^) masses are used as in- 
puts [l|. The uncertainties on these values are included 
in the systematic uncertainty on the B*2 mass. In addi- 
tion, the mass is corrected by an amount em to account 
for the DO momentum scale uncertainty. This correction 
is in proportion to the difference between the mass of the 
_B+ as measured by DO, and as listed by the PDG [l|], 
leading to an upward shift in mass eM — +0.07 MeV/c^. 
A 100% systematic uncertainty is assigned to this correc- 
tion. Taking all factors into account, the mass M{B*2) 
is measured to be: 



m{b: 



5839.6 ± 1.1 ± 0.7 MeV/c^ 



(7) 



where the first uncertainty is statistical, the second sys- 
tematic. 

Taking the detected number of B+ (20915 ± 293) and 
B*2 (125 ± 25) candidates, the production rate of i?*2 
relative to that of i?+ is calculated as follows: 



Rj = 



Br{h ^ B*2 ^ B+K-) _ N{B: 



Br{b -^ B+) 
= (1.15 ±0.23 ±0.13)%. 



N{B+)-e 

(8) 

Here e is the relative detection efficiency of B*2 events 
compared to 5+ events, i.e. it is the efficiency to select 
the additional kaon from the 5*2 decay. The value of 
this parameter is determined from simulation to be e = 
0.518 ± 0.011 (stat.), where the uncertainty results from 
the finite size of the simulation and is thus propagated 
into the measurement of Rj as a systematic uncertainty. 
Emphasis is placed on agreement between the transverse 
momentum distributions in data and in simulation, and 
a systematic uncertainty is assigned to e to account for 
any difference. 

Theoretical models predict that the 5*2 meson, exclud- 
ing phase-space factors, should decay with equal branch- 
ing ratios into B*+K and B^K. Decays into B*+ K will 
be observed as a resonance displaced to lower AM by the 
missing photon energy 45.78 ±0.35 MeV [1|. An observa- 
tion of this kind has already been made with the excited 
states of the ibd) quark system [9|] . 




M(B*K") - M(B*) - M(K") (MeV/c^) 

FIG. 2: Invariant mass difference AM = M{B+K') - 
M{B^) ~ M{K~) for exclusive B decays. The line shows 
the fit with a two-peak hypothesis, as described in the text. 
Shown separately are contributions from signal and back- 
ground. 



Since the mass difference in the decay 3*^2 ~^ B*^K 
is very small, the rate should be strongly suppressed by 
a factor proportional to {P* /P)^, where P* (P) is the 
momentum in the center-of-mass frame of the kaon in the 
decay B*^ ^ B*+K {B+K) 0. Using the B*^ mass as 
measured here, a suppression factor of 0.074 is calculated; 
therefore no detectable B*2 -^ B*^K signal is expected 
in the AM distribution with the current statistics. 

To test for the presence of a Bsi signal in the data, a 
two-peak hypothesis is used to fit the AM distribution. 
The Bsi peak is assigned a physical width of MeV/c^, 
and parameterized by a double-Gaussian function repre- 
senting the experimental detector resolution. The reso- 
lution parameters are fixed from a separate simulation of 
Bsi -^ B*^K~ events. In this case, the widths <Ji^2{Bsi) 
of the narrow and wide Gaussians are determined to be 
1.1 and 2.2 MeV/c^ respectively, and the normalisation of 
the narrow Gaussian is 3.6 times that of the wide Gaus- 
sian. Again, the widths of the Gaussians are increased 
by 10% to correct for underestimation in simulation. 

The resulting fit is shown in Fig. [21 giving the following 
parameters for the Bsi signal: 

AM{Bsi) = M{B+K-) ^ M{B+) ^ M{K-) 

= M{Bsi) - M{B*+) ~ M{K-) 

^ 11.5 ±1.4 (stat.) MeV/c^ 

N = 25 ±10 (stat.) events. (9) 

Without the Bsi signal contribution, the log likelihood of 
the fit decreases by 2.7, which implies that this structure 
is observed with a statistical significance of less than 3tT. 
Hence with the current data, the existence of a Bsi state 
can be neither confirmed nor excluded. 

For the i?*2 mass fit, the influences of different sources 
of systematic uncertainty are estimated by examining the 



changes in the fit parameters under a number of vari- 
ations. The parameters describing the background are 
allowed to vary in the fit and their uncertainties are in- 
cluded in our results. The effect of binning is tested by 
varying the bin width and position. In addition, the fit 
is made without the 10% mass resolution correction. To 
check the effect of fixing the physical width F of 3*2 at 
1.0 MeV/c^, the fit is repeated with different widths in 
the range 0-2 MeV/c^. The uncertainty in the abso- 
lute momentum scale, which results in a small shift of all 
measured masses, is assigned a 100% systematic uncer- 
tainty. Finally, the uncertainties on the PDG masses of 
-B"*" and K~ [IJ are propagated into the systematic uncer- 
tainty on the -6*2 mass. The summary of all systematic 
uncertainties in the B*2 mass fit is given in Table HI 



TABLE I: Systematic uncertainties of the B*2 parameters de- 
termined from the AM fit and from the conversion into the 
mass A'I{Bg2)- The rows show the various sources of system- 
atic uncertainty as described in the text. The columns show 
the resulting uncertainties for the two free signal parameters 
as described in Eq. Q. 



Source 



5M(B;2) (MeV/c^ 



5N 



in 


widths/positions 


0.3 


7 


alu 


eof r 


0.3 


5 


DC 


t mass uncertainties 


0.5 





Ion 


lentum scale 


0.1 





.eso 


lution uncertainty 


0.1 


3 



Total 



0.7 



The measurement of the relative production rate Rj 
uses the kaon detection efficiency predicted in simula- 
tion, as well as the numbers of 3*2 and B^ events. The 
systematic uncertainty on the number of B~^ events, de- 
scribed in Ref. [9[, is ±200 events. The systematic un- 
certainty on the number of 3*2 events is ±9 events (see 
Table U. 

The uncertainty of the impact parameter resolution 



14j . It can 



in the simulation is estimated to be ~10% 
influence the measurement of the selection efficiency of 
the kaon from the 3*2 decay. To test for the effect of 
such an uncertainty, the efficiency is recalculated with 
the kaon impact parameter requirement varied by ±10%. 
The resulting variation in efficiency is ±0.022. 

The track reconstruction efficiency for particles with 
low transverse momentum is measured in Ref. 15|, and 



good agreement between data and simulation is found. 
This comparison is valid within the uncertainties of 
branching fractions of different 3 semileptonic decays, 
which is about 7%. This uncertainty translates to an 
efficiency variation of ±0.036. An additional systematic 
effect, associated with the difference in the momentum 
distributions of selected particles in data and in simula- 
tion, is taken into account. This yields an uncertainty in 
the efficiency of ±0.002. 

Combining all these effects in quadrature, the total sys- 



TABLE II: Systematic uncertainties in the _B*2 production rate measurement. The rows show the various sources of systematic 
uncertaintys as described in the text. The columns show the effect of these sources on the three parameters used in the Rj 
measurement, and on the production rate itself. 



N{B*2) uncertainty 
N{B^) uncertainty 
Reweighting correction 
Impact parameter resolution 
Track reconstruction efficiency 
Statistical effects from simulation 



S[N{Bt2)] 



miB+)] 



S{e) 



S{Rj)i%) 



9 



200 



0.002 
0.022 
0.036 
0.011 



0.08 
0.01 
0.00 
0.05 
0.08 
0.02 



Total 



200 



0.044 



0.13 



tematic uncertainty on the efficiency e is 0.042. Both this 
and the statistical uncertainty 0.011 on e must be propa- 
gated into the production rate measurement. The effects 
of contributions from the efficiency, and the number of 
detected B'^ and B*2 candidates, are shown in Table HIl 

Different consistency checks of the observed signal are 
performed. Variation of the selection requirements on 
Pt{K), M{B+) and Sk, from B;^ -^ B+R- decays, 
does not produce any significant change in the results. 
Events with positively and negatively charged kaons are 
analyzed separately, and consistent results are obtained. 
A complementary sample of events containing a kaon not 
compatible with the primary vertex is selected, and no 
significant B*2 signal is observed. Events with wrong 
charge combinations {B~^K^ and B^K^) also show a 
signal consistent with zero. 

In conclusion, the B*2 state is observed in decays 
to B^K'' with a statistical significance of more than 
5(T. The measured mass is 5839.6 ± l.l(stat.)±0.7(syst.) 
MeV/c^. The i?*2 relative production rate with respect 
to the B+ meson is [1.15 ± 0.23(stat.)±0.13(syst.)]%. 
Searching for a i?si signal gives inconclusive results with 
the currently available data set, which is expected to in- 
crease by a factor of five in the next few years. 
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